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Resumen 

La dinámica de fluidos tiene como objetivo comprender el movimiento 
de líquidos y gases mediante funciones que describen la distribución de 
velocidades. En el presente estudio se investigan las características de 
las corrientes secundarias en canales trapezoidales con ángulos laterales 
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de 45, 60 y 75% sometidas a la influencia de cinco descargas diferentes 
con pruebas experimentales; los resultados se comparan con los de 
Flow3D. Los resultados obtenidos de las mediciones experimentales y la 
comparación de modelos numéricos de la velocidad de las corrientes 
secundarias en diferentes números de Froude demostraron que existe 
una relación opuesta entre la velocidad secundaria en la dirección 
perpendicular al eje de flujos (V,) y la velocidad en una dirección 
perpendicular al nivel de flujo (V,) en canales trapezoidales. Además, a 
un ángulo de 45%, ha habido una pérdida de energía notable durante el 
salto hidráulico. La relación del aumento en la velocidad de las 
corrientes secundarias en la dirección X en Froude número 10 de ángulo 
de 45% es mayor que la de otras dos secciones, que es igual a 71%, en 
comparación con el ángulo de 75% en modelos numéricos; en Froude 
número 9 fue de 91% durante las pruebas experimentales. Luego, la 
velocidad de las corrientes secundarias en la dirección Z para Froude 
número 2 de 75% es mayor que la de otras dos secciones, que es 88%, 
en comparación con el ángulo de 45% en modelos numéricos y en Froude 
número 1.5 fue igual a 74.5% para pruebas experimentales. 


Palabras clave: corrientes secundarias, canal trapezoidal, salto 
hidráulico, tensión de cizallamiento límite. 


Abstract 


Fluid dynamics aims at understanding the movement of liquids and 
gases by functions that describe the distribution of velocities. In the 
present study, the characteristics of secondary currents in trapezoidal 
channels with side angles of 45, 60 and 75%, subjected to the influence 
of five different discharges with experimental tests, are investigated and 
the results are compared with that of Flow3D. The results obtained from 
experimental measurements and numerical models comparison of the 
rate of secondary currents in different Froude numbers demonstrated 
that there exists an opposite relationship between the secondary 
velocity in the direction of perpendicular to the axis of flows (V,) and 
velocity in a direction perpendicular to the flow level (V,) in trapezoidal 
channels. Moreover, at 45% angle, there has been a remarkable energy 
loss during hydraulic jump. The ratio of the increase in the secondary 
currents velocity in X direction in Froude number 10 of 45% angle is 
higher than that of two other sections, which is equal to 71%, compared 
to 75% angle in numerical models and in Froude number 9 it was 91% 
during experimental tests. Then, the secondary currents velocity in Z 
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direction for Froude number 2 of 75% angle is higher than that of two 
other sections, which is 88%, compared to the 45% angle in numerical 
models and in Froude number 1.5 it was equal to 74.5% for 
experimental tests. 


Keywords: Secondary currents, trapezoidal channel, hydraulic jump, 
boundary shear stress. 


Recibido: 31/07/2018 
Aceptado: 13/12/2018 


Introducción 


Uno de los casos más frecuentes de flujo rápidamente variado es el 
fenómeno del salto hidráulico. La aplicación más importante del salto 
hidráulico es en la disipación de energía debajo de vertederos, 
vertederos, puertas, etcétera. Basado en observaciones en 3D durante 
pruebas experimentales, los fenómenos de flujo difieren de manera 
significativa del salto hidráulico clásico en canales trapezoidales y 
rectangulares. El caso más estudiado del salto hidráulico es en canales 
rectangulares horizontales, a través de discusión detallada métrica. El 
dibujo esquemático de un salto hidráulico típico en un canal horizontal 
se muestra en la Figura 1. El salto hidráulico y su estabilidad en el canal 
trapezoidal son especialmente importantes. Académicos como 
Wanoschek y hHager (1989), y .Muhsun (2012) analizaron las 
propiedades de los canales trapezoidales a través de pruebas 
experimentales. Hager (1992) y Chanson (2007) estudiaron el salto 
hidráulico unular y describieron sus características, donde se calcularon 
los valores del número de Froude en el que el salto no es más que un 
movimiento irregular, con independencia del efecto del ancho del canal. 
Afzal y Bushra (2002) estudiaron la estructura del salto hidráulico 
turbulento en el canal trapezoidal. 
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Figura 1. Estructura de flujo de salto hidráulico. 


El perfil de la superficie se aproxima como una solución universal 
limitante siempre que las variables estén apropiadamente no 
dimensionadas. Gupta, Mehta y Dwivedi (2013) investigaron la longitud 
relativa y la pérdida relativa de energía del salto hidráulico libre en 
canales prismáticos. Ohtsu (1976) clasificó los saltos hidráulicos en 
canales trapezoidales en cuatro tipos. El criterio implicaba la simetría de 
la corriente principal. Ohtsu (1976) recomendó que sólo se usen saltos 
hidráulicos de los tipos 1 y 2 ("el salto es casi simétrico") como 
disipador de energía. Kim, Choi, Park y Byeon (2015) han investigado 
las propiedades del salto hidráulico y su disipación de energía en aguas 
abajo de la compuerta. Javan y Eghbalzadeh (2013) han estudiado salto 
hidráulico sumergido con modelado turbulento K-e. Nezu et al. (1985), 
después de estudiar el efecto de la rugosidad en los resultados 
experimentales en un canal rectangular, incluso cuando las condiciones 
de rugosidad límite cambian, mostraron que la estructura principal de 
las corrientes secundarias no se modifica de forma drástica. Yaaghoubi y 
Givechi (2007) investigaron el efecto de las corrientes secundarias en la 
velocidad promedio de las secciones transversales en canales 
rectangulares. Después de estudiar el efecto de la rugosidad en los 
resultados experimentales en un canal rectangular, incluso cuando 
cambian las condiciones de la rugosidad de los límites, la estructura 
principal de las corrientes secundarias no cambia dramáticamente. 
Nezu, Nakagawa y Tominaga (1985), y Qin, Shao y Zhou (2016) 
hicieron una comparación entre dos lunares de corrección de corriente 
secundaria diferentes para el flujo de profundidad media de los canales 
serpenteantes. Otros investigadores, como Dasgupta y Tomar (2015), 
también investigaron la distribución de las tensiones de cizallamiento 
límite en canales trapezoidales prismáticos. Shokrian-Hajibehzad y 
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Shafai-Bejestan (2017) predijeron la relación de profundidad secuencial 
en el salto hidráulico en el río y encontraron que la relación de 
profundidad secuencial en lechos rocosos ásperos es considerablemente 
más pequeña que en una cama de concreto lisa. Se presentan dos 
nuevas ecuaciones de predicción para pendientes de lechos suaves y 
ásperos; los resultados mostraron que las predicciones concuerdan de 
forma razonable con los datos experimentales. 


Jalili-Ghazizadeh, Attari y Farhadi-Rad (2017) estudiaron el salto 
hidráulico junto con presas. En el modelo se considera el momento de 
salida lateral; el modelo predice la ubicación del salto hidráulico. El 
presente estudio busca llenar el vacío de investigación relacionado con 
la investigación de células de corriente secundaria en canales 
trapezoidales y su efecto en el fenómeno del salto hidráulico, que debe 
considerarse como la novedad de la investigación. 


El presente documento persigue dos objetivos principales: 1) las 
características y propiedades de las corrientes secundarias en el salto 
hidráulico para canales trapezoidales que incorporan tres diversos 
ángulos de inclinación lateral (ángulos laterales), que utilizan enfoques 
experimentales y numéricos, y 2) variación en la pérdida de energía en 
el salto hidráulico en relación con los cambios en los números de Froude 
para los análisis tanto experimentales como numéricos, así como la 
comparación con los hallazgos obtenidos por Hager (1992). El enfoque 
principal del estudio se concentra en diferentes secciones trapezoidales 
sometidas a varios flujos de descarga. 


Células de corriente secundaria 


En general, las corrientes secundarias son un subconjunto de flujos en 
las secciones transversales o alrededor del eje perpendicular a la sección 
transversal del flujo. Las corrientes en las curvas se llaman corrientes 
secundarias fuertes, pero las que se encuentran en los canales 
prismáticos del plan se denominan corrientes secundarias débiles. Las 
corrientes secundarias y los canales abiertos son significativamente 
importantes en la ingeniería hidráulica debido al inmenso efecto de 


Tecnología y ciencias del agua, 10(3), 190-218. DO!: 10.24850/j-tyca-2019-03-08 


194 


2019, Instituto Mexicano de Tecnología del Agua 
Open Access bajo la licencia CC BY-NC-SA 4.0 
Tecnología y (https: //creativecommons.org/licenses/by-nc- 


CienciaszAgua sa/4.0/) 


estas corrientes en el campo de flujo, el impacto de la regla de fricción y 
la formación de formas tridimensionales en los lechos de los canales. 


Ecuaciones de corrientes secundarias 


El propósito principal de hacer corrientes secundarias en un flujo 
longitudinal uniforme es la homogeneidad de la disparidad de tensión de 
Reynolds y, en esencia, de sus componentes verticales. Prandtl (1904) 
dictó que los dos mecanismos principales conducen a flujos de rotación 
en la dirección del flujo, es decir: 1) desviación de las tensiones de 
cizallamiento principales por un gradiente de presión lateral o fuerza 
física; 2) inconsistencias y turbulencia fluctuante en una placa 
perpendicular a la del eje de rotación. La ecuación de rotación 
longitudinal en el flujo turbulento completamente desarrollado está 
subrogada (Azhdari-Moghaddam, Tajnasaj €: Givech, 2013): 


yy? (0 -w)+ 

Oy Oz  0OyOz (1) 
o 2 o e 

A AT — 

o id Gr al 

Ov 0w 

ES (2) 
Oz  0y 


Donde W y V se refieren a componentes de velocidad en dirección 
horizontal (Z), y componentes de velocidad en dirección vertical (Y), 
respectivamente. 


De hecho, estos son los mismos componentes de la corriente 
secundaria, donde Y es la viscosidad cinemática del fluido, denota el 
esfuerzo cortante de Reynolds y se refiere al estrés normal de Reynolds. 
Eso está indicado por la función de flujo: 
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La Figura 2 ilustra la diferencia de tensiones verticales en la estructura 
principal de los flujos secundarios. 


172 172 
Figura 2. Líneas no dimensionales (w"-v )) 27.4. desde la derecha en 
U*máx 


el orden de canal cerrado, canal trapezoidal y canal rectangular. 


En el presente trabajo, las variables Z e Y son del orden de X para la 
velocidad perpendicular al eje de flujo y Z para la velocidad 
perpendicular al lecho del canal. La siguiente ecuación muestra el 
fenómeno del salto hidráulico en el estudio actual, según se investigó, 
en función de la definición de los flujos de rotación, de modo que la 
curvatura funciona en flujo opuesto a cero: 


io jo k 

SA 0 0 0 a 

VxU= | a eijkOjvk + 0 (4) 
Ux Uy Uz 


Donde V representa el operador nebula, calculado como sigue en el 
sistema cartesiano. 


U es componente de velocidad: 
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2..0 . 0 0 

1 (5) 


Modelado de flujo 3D 


La resolución del problema de CFD (dinámica de fluidos computacional), 
por lo general consta de cuatro componentes principales: geometría y 
generación de cuadrícula, configuración del modelo físico, solución y el 
proceso posterior de los datos computados. El problema del conjunto se 
calcula para generar la geometría y la cuadrícula, y los datos se 
adquieren de una manera bien conocida. El modelado de un fenómeno 
implica problemas complejos que pueden resolverse desarrollando un 
modelo simple. Por lo tanto, aunque un modelo ideal debería introducir 
la cantidad mínima de complejidad en las ecuaciones de modelado, es 
importante ser lo más preciso posible. Para llevar a cabo un modelo 
numérico de las corrientes secundarias, inicialmente es necesario 
desarrollar un modelo numérico de las secciones indiferentes de la 
velocidad media para obtener los límites de un cambio de velocidad de 
forma adecuada. 


Según Bayon, Valero, Garcia-Bartual y López-Jiménez (2016), se suele 
suponer que el valor límite para las cantidades KT y eT determinan el 
equilibrio local entre los procesos de producción y caída de la tensión de 
corte turbulenta y una regla de la velocidad de la pared perfil. En esta 
investigación, el modelo RNG (grupo de renormalización) se usó debido 
a sus métodos estadísticos para derivar ecuaciones promedio para 
cantidades de perturbación, como la energía cinética y su tasa de 
pérdida, lo que lo convierte en el mejor modelo de turbulencia para 
simular un salto hidráulico. El software Flow-3D se emplea para obtener 
una correlación precisa entre la solución y convergencia de respuestas, 
que se basan en el flujo espectral geométrico y las características del 
flujo (Microsoft Corporation, 2015). 


Correlación de la señal ADV 
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Se pueden encontrar explicaciones extensas sobre los principios de 
operación de ADV en el Manual Operativo de Nortek (Nortek-AS, 2013). 
Las velocidades se estiman utilizando una técnica denominada 
procesamiento Doppler coherente de pulso, donde el ADV (velocímetro 
doppler acústico) mide el cambio en la fase de la señal de retorno a 
partir de dos pulsos acústicos sucesivos. El retorno acústico no es un 
reflejo de un solo objetivo sino una superposición de los reflejos de 
muchas partículas individuales contenidas en el volumen de muestreo. 
Se logra una coherencia de fase si todas las partículas en el volumen de 
muestreo mantienen sus posiciones relativas entre sí, en cuyo caso la 
fuerza y las fases relativas de las reflexiones individuales no cambiarían 
de un pulso al siguiente. Sin embargo, el ruido doppler es una parte 
inherente de los sistemas de retrodispersión de volumen basados en 
doppler, y también existen otras fuentes de ruido. Por lo tanto, la señal 
de retorno del segundo pulso no es una reproducción de cambio de fase 
del primer pulso, sino que contiene una cierta cantidad de ruido, que se 
puede agregar a la parte coherente de la señal de retorno (Nortek-AS, 
2013): 


Donde Si 'y S2' son las señales de retorno de dos pulsos adyacentes; S; 
y S2 son las partes coherentes de la señal, y N; y N, representan el 
ruido aleatorio. En la práctica, el ADV calcula el cambio de fase entre las 
señales de retorno de dos pulsos sucesivos utilizando una función de 
autocorrelación compleja, que tiene una fase y una magnitud. La fase es 
proporcional a la distancia que se recorre en el tiempo entre los pulsos 
por dispersión y, por lo tanto, a la velocidad. La magnitud es una 
medida de qué tan similares son los ecos de dos señales de retorno, y la 
magnitud normalizada produce el coeficiente de correlación. Si los ecos 
de dos pulsos vecinos fueran idénticos y sólo con un cambio de fase, la 
correlación sería perfecta e igual a 100%. A medida que los ecos se 
vuelven más disímiles, la correlación disminuye. Las correlaciones bajas 
se asocian con flujos altamente turbulentos, burbujas de aire en el flujo, 
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valores bajos de SNR (relación señal/ruido), gradientes de gran 
velocidad dentro del volumen de muestreo, la presencia de partículas 
individuales grandes o interferencia desde el límite. 


En flujos turbulentos, como el salto hidráulico, es razonable suponer que 
las partículas en el volumen de la muestra no mantendrán sus 
posiciones relativas entre sí, por lo que se agregará el ruido doppler de 
la señal de retorno. Esto descorrelatará la señal y, en este caso, las 
correlaciones de señal baja no significarán de forma necesaria que los 
datos sean indeseables. Además, si los remolinos de turbulencia son del 
mismo orden que el volumen de la muestra o más pequeños, la 
correlación cambiará aún más, lo que conducirá a un aumento del ruido 
en las mediciones de velocidad. El coeficiente de correlación para flujos 
turbulentos se tratará con más detalle en la siguiente sección. 


Estudio experimental 


Para crear un salto hidráulico a lo largo del canal se han probado varios 
modos hidráulicos en el modelo de laboratorio. Estudios anteriores 
examinaron los saltos clásicos hechos para secciones transversales 
rectangulares. La configuración experimental y las propiedades del canal 
se presentan en la Tabla 1. 


Tablal.Características geométricas del canal. 


Tipo de Long | Anchura Alto de Válvula, — Ubicación 
canal itud (m) pared apertura de la 
(m) lateral(m) (m) válvula(m) 
m=1.0 
m = 0.58 5 0.2 1 0.035 0.2 
m=0.26 
Donde m es la pendiente lateral del canal y se representa 


esquemáticamente en la Figura 3. El valor numérico de la pendiente 
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lateral para los ángulos 45, 60 y 75% en esta investigación fue igual a 
1.0, 0.58 y 0.26, respectivamente. 


Figura 3. Pendiente lateral en media sección transversal. 


El volumen del tanque de reserva de inicio del canal es de 3.2 m”. La 
Figura 4 y la Figura 5 muestran las características de geometría del 
canal, y la posición del ADV en el canal. 


Figura 4. Geometría del canal trapezoidal en tres ángulos. 
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Figura 5. Especificaciones del dispositivo ADV en el canal. 


Como se muestra en la Figura 5, la máquina ADV está montada en un 
carro, se mueve en el canal utilizando el software (Jog) y envía el 
dispositivo a la profundidad designada en la ubicación deseada. En este 
documento, las mediciones se llevaron a cabo a 25 Hz durante dos 
minutos, con un total de 3 000 impresiones para la velocidad media del 
punto (u, v, w). 


Basado en la experiencia previa, el dispositivo al inicio descansa por 
completo en flujo laminar, de modo que puede enviar la frecuencia 
dentro del flujo; luego se usó una descarga más baja para cosechar 
velocidades secundarias al comienzo de un salto hidráulico. Las 
descargas utilizadas en este estudio fueron iguales a 10, 30, 50, 70 y 90 
l/s, para tres ángulos de sección geométrica, incluidos los ángulos de 
45, 60 y 75% en el canal trapezoidal. 


El punto que importa en la configuración del software Victorina+ es que 
todos los valores en la configuración de prueba y error deben coincidir 
de manera coordinada con la obtenida para los cuatro destinatarios de 
que el coeficiente de correlación requiere ser mayor que 70% y la SNR 
debe ser mayor que 15 dB (Martin, Fisher, Millar, £ Quick, 2002). En el 
presente estudio se obtiene el valor numérico de SNR de 23 a 24 dB. 


Cabe señalar que el coeficiente de correlación es un parámetro para 
determinar la calidad de la medición de la velocidad en porcentaje. En 
cada medición, el dispositivo calcula el parámetro de correlación para 
cada receptor de audio. El valor de correlación de 100% indica, en el 
mejor de los casos, una medición de ruido bajo; una correlación cero 
señala el efecto de la reducción de ruido. En lo ideal, las mediciones 
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deben tener una correlación entre 70 y 100%. Además, la relación SNR 
es señal de ruido y la cantidad de este parámetro en las mediciones es 
indicativa de la presencia de partículas suficientes en el agua para 
dispersar el sonido. Si el agua es muy clara y suave, la señal de retorno 
a los receptores del dispositivo es débil, en comparación con el ruido 
existente, y el dispositivo no podrá medir la velocidad. Durante la 
medición de flujos turbulentos, como saltos hidráulicos, tal relación no 
debe ser inferior a 15. En la presente investigación, las bombas se 
utilizaron para determinar la descarga deseada en el canal. Con respecto 
a la Tabla 2, el estado de las bombas rpm está listo para la realización 
de diferentes descargas. Vale la pena señalar que durante la prueba se 
utilizaron tres bombas, de tal manera que todas las demás descargas se 
formaron mediante la combinación de 2 o 3 bombas, esperando los 
valores de las dos primeras descargas. 


Tabla 2. El estado de las bombas ultrasónicas. 


Fila Tipo de Ronda por minuto Q (l/s) 
bomba (RPM) 
1 Descarga 1 754 33 
constante 
2 4 783 90 
3 3 507 70 
4 Descarga 2 659 50 
5 variable 1 576 30 
6 896 10 


El dispositivo Vectrino mide la velocidad del agua según el fenómeno 
doppler. Basado en este fenómeno, si la fuente de audio se mueve a 
una velocidad v a un receptor sensible, la frecuencia de audio recibida 
es calculada por el receptor de acuerdo con la velocidad de audio del 
transmisor usando la siguiente ecuación: 


Eboppler E Fuel) (8) 


Donde Fboppler está cambiando la frecuencia del sonido recibido y Fsource 
es la frecuencia de sonido del transmisor; v es la velocidad del 
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transmisor al receptor; c es la velocidad del sonido en el fluido, que se 
supone que es 1 484 m/s en el presente estudio. 


La Figura 6 muestra cómo se coloca el dispositivo dentro de la corriente, 
que comienza a fluir a una descarga de 10 a 90 l/s; cada descarga se 
tomó durante dos minutos, y se realizó en los tres ángulos del canal. 


Figura 6. Cómo colocar el dispositivo dentro de la corriente. 


En la Figura 7, Figura 8 y Figura 9 se muestra el comienzo del salto 
hidráulico formado por la descarga de 90, 70 y 50 l/s a 75". 


Figura 7. Comienzo del salto hidráulico en 90 l/s. 
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Figura 8. Inicio del salto hidráulico en 70 l/s. 


Figura 9. Inicio del salto hidráulico en 50 l/s. 


La Tabla 3 y Tabla 4 presentan resultados experimentales y numéricos, 
respectivamente, e ilustran cinco descargas diferentes de 90, 70, 50, 30 
y 10 l/s. Estas tablas también presentan los números de Froude antes y 
después del salto, la profundidad inicial y la relación de la longitud del 
salto a la profundidad secundaria para las tres secciones geométricas, lo 
que se hace para mostrar no dimensional de la longitud del salto. 


Tabla 3. Resultados experimentales en tres secciones trapezoidales. 


m=1i.0 m = 0.58 m = 0.26 
(F-) | L; Profun- (F,) UL; Profun- (F) | L; Profundidad 
Y, didad Y, didad Y, inicial (m) 
inicial (m) inicial 
(m) 
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Antes del 
S Sto. (ar 618 120 de 12 | 5.11 
20 Después 0.06 0.073 0.08 
(/s) | del salto | y 5 0.64 197 
Antes del 
E Saló a la 7 | 53 4 6 
0 Después 0.05 0.057 0.062 
(l/s) | del salto 0.73 0.76 0.84 
Antes del 
= Salta 2 2.1 | 5.41 salsa 
an Después 0.027 0.037 0.042 
(Vs) | del salto | y 36 0.63 0.81 
Antes del 
E salto 2 4 1.86 | 4.9 1.1 | 47 
AD Después 0.024 0.03 0.038 
(Vs) | del salto | y 78 0.31 0.38 
Antes del 
E salto 12 | 3.1 1 ls 1.2 | 3.81 
en Después 0.025 0.031 0.033 
(Vs) | del salto | y ,3 0.17 0.16 
Tabla 4. Resultados numéricos en tres secciones trapezoidales. 
m=1i.0 m = 0.58 m = 0.26 
(F,) AL Profundidad | (F,) A Profundidad | (F,) AL Profundidad 
Y | inicial (m) Y | inicial (m) Y | inicial (m) 
= | Antes |9.72 10 | 6.83 0.055 9.2 |6.18 0.071 
90 | del salto 6.92 0.05 
(l/s) | Después | 0.73 0.87 0.87 
del salto 
== Antes 8.12 7.6 | 6.78 0.049 8.13 | 6.13 0.052 
70 | del salto 6.92 0.041 
(l/s) | Después | 0.93 0.98 0.99 
del salto 
= Antes 5.39 7.31 | 5.35 0.036 7 5.5 0.043 
50 | del salto 6.74 0.031 
(l/s) | Después | 0.97 0.59 0.33 
del salto 
Q= | Antes |2.97 4 (5.25 0.026 47 4 0.032 
30 | del salto 
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(l/s) | Después | 0.34 | 5.76 0.02 0.75 0.94 
del salto 
= Antes | 2.48 2.62 | 4.4 0.033 3.76 | 3.1 0.039 
10 | del salto 4.9 0.021 
(l/s) | Después | 0.2 0.93 0.5 
del salto 


La Figura 10a y la Figura 10b representan la dispersión de SNR y 
correlación, respectivamente, en corrientes secundarias en la dirección 
perpendicular al eje de flujo (V,) en relación con el número de 
impresiones del dispositivo en el salto hidráulico. La Figura 11a y la 
Figura 11b muestran la dispersión de SNR y la correlación en corrientes 
secundarias en la dirección perpendicular al eje de flujo (V,) en relación 
con el número de impresiones del dispositivo en el salto hidráulico. Con 
respecto a la Figura 10 y la Figura 11, la distribución de la correlación a 
las velocidades de las corrientes secundarias en la dirección 
perpendicular al eje (V,) del flujo fue mayor que las velocidades 
secundarias en la dirección de la perpendicular al lecho del canal (V,). 
En altas descargas, debido a la mayor velocidad de flujo, la correlación 
se mitiga en la dirección perpendicular al eje de flujo (V,). 


25 , 
24 - 
2 23, $ 
o % + + 


22 

y? 
21 + r 
0 1000 2000 3000 4000 
Pasos de tiempo basados en la frecuencia de ADV 


Correlación 
mn 
o 


0 1000 2000 3000 4000 
Pasos de tiempo basados en la frecuencia de ADV 


Figura 10. (a) El proceso de adquisición de datos cambia en función de los pasos de 
tiempo de la frecuencia de ADV en relación con la SNR en las corrientes secundarias en 
la dirección V,; (b) el proceso de adquisición de datos cambia en función de los pasos 
de tiempo de la frecuencia de ADV en relación con la correlación de la señal en 
corrientes secundarias en la dirección V,. 
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4 + 
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Figura 11. (a) El proceso de adquisición de datos cambia en función de los pasos de 
tiempo de la frecuencia ADV en relación con el SNR en corrientes secundarias en la 
dirección V,; (b) el proceso de adquisición de datos cambia según los pasos de tiempo 
de la frecuencia de ADV en relación con la correlación de la señal en corrientes 
secundarias en la dirección V,. 


La tendencia de la distribución de la velocidad en la dirección horizontal 
(9) fue mayor en los canales con menos pendientes de la pared lateral, 
como las que tienen un ángulo de 45". La velocidad del flujo también es 
mayor en la dirección horizontal (x), por lo tanto conduce al desarrollo 
de la tensión de corte en los canales y finalmente reduce la energía del 
salto hidráulico en la geometría del canal, que en efecto tiene un mayor 
efecto en descargas más altas. En línea con los objetivos de esta 
investigación y para estudiar las celdas de la corriente secundaria en la 
dirección horizontal (x), se analizaron dos altas descargas: 90 y 70 l/s. 
Las descargas por debajo de estos valores no afectan de forma 
dramática la pérdida de energía de salto en las celdas de corriente 
secundarias en la dirección horizontal (x). Con los aumentos en la 
velocidad de flujo a lo largo del canal (y), la velocidad de las corrientes 
secundarias en la dirección x aumenta más en un ángulo de 45%, en 
comparación con los ángulos de 60 y 75%. 


En contraste, la distribución del vector de velocidad de las corrientes 
secundarias perpendiculares a la superficie de flujo (z) en la descarga 
más débil es mayor que la de las paredes laterales, debido a la baja 
velocidad del flujo y al predominio de la velocidad de distribución de la 
gravedad en el piso del canal. Este es el caso en los valores de velocidad 
de los flujos secundarios en cada una de las tres secciones que se 
muestran en la Tabla 5 asociadas con las pruebas experimentales y en 
la Tabla 6 relacionadas con el modelado numérico. Como se muestra en 
la Tabla 5 y Tabla 6, los valores de velocidad están en la dirección 
horizontal (x) para altas descargas, y los valores para los ángulos de 45 
y 60% están en la dirección vertical. Además, la velocidad de las 
corrientes secundarias en la dirección vertical (z) para descargas bajas 
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en ángulos de 60 y 75% fue mayor que la de las corrientes secundarias 
en la dirección horizontal (x). 


Tabla 5. Velocidades de flujo de corrientes secundarias experimentales 
en x y z directas. 


Fila Q (l/s) m=1i.0 m = 0.58 m = 0.26 
1 90 V, = 0.36 m/s V, = 0.31 m/s V, = 0.26 m/s 
V,= 0.22 m/s V,= 0.25 m/s V,= 0.27 m/s 
2 70 V, = 0.33 m/s V, = 0.18 m/s V, = 0.15 m/s 
V,= 0.26 m/s V,= 0.29 m/s V, = 0.49 m/s 
3 50 V,= 0.21 m/s V,= 0.16 m/s V,= 0.11 m/s 
V,= 0.39 m/s V, = 0.5 m/s V, = 0.6 m/s 
4 30 V,= 0.14 m/s V,= 0.1 m/s V,= 0.07 m/s 
V,= 0.41 m/s V,= 0.52 m/s V,= 0.81 m/s 
5 10 V,= 0.11 m/s V,= 0.05 m/s V,= 0.01 m/s 
V, = 0.49 m/s V,= 0.61 m/s V,= 0.92 m/s 


Tabla 6. Velocidades de flujo de corrientes secundarias numéricas en x 


y z directas. 
Fila | Q (l/s) m=1.0 m = 0.58 m = 0.26 
1 90 V,= 0.39 m/s V, = 0.35 m/s V,= 0.29 m/s 
V,= 0.19 m/s V,= 0.23 m/s V,= 0.25 m/s 
2 70 V,= 0.36 m/s V,= 0.21 m/s V,= 0.17 m/s 
V,= 0.24 m/s V,= 0.27 m/s V,= 0.46 m/s 
3 50 V,= 0.25 m/s V, = 0.18 m/s V,= 0.15 m/s 
V, = 0.35 m/s V, = 0.46 m/s V,= 0.57 m/s 
4 30 V,= 0.16 m/s V, = 0.13 m/s V,= 0.1 m/s 
V, = 0.36 m/s V,= 0.47 m/s V,= 0.76 m/s 
5 10 V, = 0.14 m/s V, = 0.11 m/s V, = 0.05 m/s 
V, = 0.41 m/s V, = 0.55 m/s V, = 0.87 m/s 
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La Figura 12 ilustra la variación en la velocidad de la corriente 
secundaria en la dirección del eje de flujo vertical (V,), que se 
representa con cambios en los números de Froude. La velocidad de las 
celdas de la corriente secundaria para la sección transversal del ángulo 
de 45% en la dirección x a diferentes tasas de número de Froude (Figura 
12) fue mayor que la velocidad de otras secciones en las pruebas 
experimentales y en el modelado numérico. Se puede concluir que la 
tensión de corte creada en las paredes laterales de esta sección es 
mayor que la tensión de corte creada en las otras secciones, lo que 
conduce el flujo hacia las paredes laterales. La Figura 12 ilustra la 
variación en la velocidad de las corrientes secundarias como se muestra 
en la dirección vertical (V,) con respecto al número de Froude. Los 
resultados muestran que la tendencia es opuesta a la dirección 
horizontal (V,). La Figura 13 ilustra los cambios en la velocidad de las 
corrientes secundarias, como se muestra en la dirección vertical (V¿) con 
respecto al número de Froude. Los resultados muestran que la 
tendencia es opuesta a la dirección horizontal (V,). 


45 
40 
al 35 
S 30 
Z 25 
E 
2 a A m=1.0 Flow3D 
> 15 * m=1.0 Experimental 
u + m=0.58 Flow3D 
10 - m=0.58 Experimental 
0 m-0.26 Flow3D 
5 + m=0.26 Experimental 
0 
0 2 4 6 8 10 12 14 16 
Froude number (Fr,) 


Figura 12. Comparación de la velocidad de las corrientes secundarias 
en horizontal directo (V,) con la función de F, para varias m en Flow 3D 
y pruebas experimentales. 
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Figura 13. Comparación de la velocidad de la corriente secundaria en 
directa vertical (V¿) con la función de F; para varias m en el caudal 
Flow3D y pruebas experimentales. 

De acuerdo con los resultados de la Figura 12, las ecuaciones que rigen 
el proceso de variación de la velocidad de las corrientes secundarias en 
la dirección horizontal (V,) frente al número de Froude para cada una de 

las tres secciones siguen las siguientes relaciones: 

Por 

m=10 Vb rmay) Xx 100 =13.61n(Fr1) +7.1 (9) 

m = 058 (V*/ 7.0) x 100 = 9.6 In(Fr1) + 6.6 (10) 

m = 0.26 (V%/7,0y) x 100 = 9.4 In(Fr1) + 3.4 (11) 

Además, de acuerdo con los resultados de la Figura 13, las ecuaciones 
que rigen el proceso de variación de la velocidad de las corrientes 
Tecnología y ciencias del agua, 10(3), 190-218. DO!: 10.24850/j-tyca-2019-03-08 210 


2019, Instituto Mexicano de Tecnología del Agua 
Open Access bajo la licencia CC BY-NC-SA 4.0 
Tecnología y (https: //creativecommons.org/licenses/by-nc- 


CienciaszAgua sa/4.0/) 


secundarias en la dirección vertical (V,) frente al número de Froude para 
cada una de las tres secciones siguen las siguientes relaciones: 


Por 

m=10 (2/70) X 100 =-11In(Pr1) + 46.4 (12) 
m=058 Vray) Xx 100 =-—14In(Fr1) + 57.7 (13) 
m =0.26 (V2/7a0)x 100 =-24.8 In(Fr1) + 86.8 (14) 


La velocidad aumenta en direcciones horizontales cuando la sección 
transversal cambia de ángulos de 75 a 45%, en números altos de 
Froude. También es posible aumentar la velocidad vertical (V,) 
cambiando la sección transversal de ángulos de 45 a 75%, en números 
bajos de Froude. Con base en los datos presentados por la Figura 12 y 
Figura 13, y la Tabla 5 y Tabla 6, se pueden obtener tasas de 
crecimiento aproximadas del 33 y 28.9% en una dirección horizontal 
(V,) de ángulos de 75 a 60% en pruebas experimentales, y en la misma 
situación en el modelado numérico, respectivamente; mientras que se 
puede obtener aproximadamente 44% en ángulos de 60 a 45“ en las 
pruebas experimentales, y 32.6% en la misma situación en el modelado 
numérico. Las corrientes de velocidad secundarias en direcciones 
verticales (V,) aumentaron en un 22.5% desde ángulos de 45 a 60% 
durante las pruebas experimentales; mientras que las corrientes 
aumentaron en un 28% durante el modelado numérico, y subieron en 
un 42% desde 60 a 75% ángulos durante pruebas experimentales. 
Además, las corrientes se incrementaron en un 47% durante el 
modelado numérico. Los resultados de los datos de la Tabla 5 y Tabla 6 
se compararon utilizando el error cuadrático medio (RMSE) y el 
coeficiente de determinación (DC), que se calcula utilizando las 
siguientes dos relaciones. Los resultados se dan en la Tabla 7. Se puede 
deducir que el efecto máximo de la velocidad de la corriente secundaria 
emergió en la dirección horizontal (V,) a 45%, y el efecto máximo de la 
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velocidad de la corriente secundaria fue en la dirección vertical (V,) en 
ángulo de 75%: 


RMSE = NE E) (15) 
_é€ y2 
DO 1 a si 


Tabla 7. Resultados de RMSE y DC para tres secciones de geometría. 


v, v, 
Tipo de DC RMSE DC RMSE 
canal 
trapezoidal 
m=1.0 0.77 0.098 | 0.91 0.067 
m = 0.58 0.89 0.084 | 0.82 0.08 
m= 0.26 0.95 0.08 0.85 0.07 


El rango de la DC está entre O y 1; a medida que su valor se aproxima a 
1, la situación se vuelve deseable. El valor de la RMSE está más cerca 
del número O tal como está. Por lo tanto, los números 0.91 y 0.067 para 
la velocidad de las corrientes secundarias en la dirección horizontal (V,), 
y 0.95 y 0.08 para la velocidad de las corrientes secundarias en la 
dirección vertical (V,) son los mejores valores posibles para alcanzar los 
resultados deseados. La Figura 14 muestra la relación de profundidades 
secuenciales Y2/Y1 como funciones de F; para varios (m) tanto en el 
modelado numérico como en las pruebas experimentales que se 
compararon con los resultados experimentales obtenidos por Hager 
(1992). La Figura 15 ilustra la comparación de la longitud de salto 
hidráulico Lj/Y1 en función de Fi para varios (m) en el trabajo actual y 
se comparó con los resultados experimentales de Hager (1992). La 
Figura 16 muestra las variaciones de la pérdida relativa de energía en 
un volcán hidráulico. En la Figura 16, la pérdida de energía en los 
modelos realizados se muestra debajo de los cambios en el número de 
Froude. Es evidente que tanto en el modelado numérico como en las 
pruebas experimentales la pérdida de energía en las secciones 
trapezoidales con un ángulo lateral de 45% fue significativa, en 
comparación con las otras dos secciones con el aumento del número de 
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Froude. Los resultados obtenidos se compararon con los de Hager 
(1992) para una mejor validación de los datos para este gráfico de 
pérdida de energía en cada una de las tres secciones de este estudio. 
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Figura 14. Secuencia de profundidad de secuencia (Y2/Y1) como 
funciones de F; para varias m en Flow 3D y pruebas experimentales. 
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Figura 15. La relación entre la longitud del salto hidráulico y la 
profundidad inicial (Lj/Y1) como una función de F, para varios m en Flow 
3D y pruebas experimentales. 
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Froude Number(F,) 


Figura 16. Comparación de la pérdida de energía relativa en el salto 
hidráulico en función de FF; para varios m en modelos Flow3-D y pruebas 
experimentales. 


Con respecto a la Figura 15, la longitud del salto hidráulico en ambos 
modelos numéricos y pruebas experimentales a un ángulo de 45% se ha 
incrementado con el aumento en el número de aterrizaje, en 
comparación con las otras dos secciones. 


Como se muestra en la Figura 16, la pérdida de energía en el salto 
hidráulico en un ángulo de sección transversal de 45% durante las 
pruebas experimentales mostró un aumento de 4%, en comparación con 
los resultados de Hager. Además, durante el modelado numérico en un 
ángulo de 45%, se observó un aumento de 9% en comparación con los 
resultados de Hager. En términos generales, estos resultados indican la 
verificación de las pruebas numéricas y experimentales. De acuerdo con 
la comparación realizada entre los resultados del modelado numérico y 
las pruebas experimentales, existe un acuerdo deseable entre los 
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resultados de estos dos modelos en las tres secciones geométricas. De 
modo que, a un ángulo de 45% en los modelos numéricos hay un 
aumento de 5%, en comparación con las pruebas experimentales, y a 
un ángulo de 60% en los modelos numéricos hay un incremento de 11% 
en comparación con las pruebas experimentales, y también a un ángulo 
de 75% en los modelos numéricos hay un aumento de 14%, en 
comparación con las pruebas experimentales. Así que el mayor 
incremento puede atribuirse a la sección de 45". 


Conclusiones 


1. Las condiciones de contorno utilizadas en este estudio son muy 
efectivas para determinar los límites de los cálculos y las sensibilidades. 


2. La velocidad de las corrientes secundarias perpendiculares al eje de 
flujo (V,) en Froude número 10 de un ángulo de 45% es mayor que las 
otras dos secciones, que es de 71%, en comparación con el ángulo de 
75% en modelos numéricos y en Froude número 9. Fue de 91% en 
pruebas experimentales. 


3. La velocidad de las corrientes secundarias en la dirección vertical de 
la superficie de flujo (V,) en el número 2 de Froude de 75% es mayor que 
las otras dos secciones, que es de 88%, en comparación con el ángulo 
de 45% en modelos numéricos y en número de Froude 1.5 fue de 74.5% 
en pruebas experimentales. 


4. Los resultados demostraron la existencia de una relación inversa 
entre la velocidad de la corriente secundaria en la dirección del eje 
vertical del flujo (W,) y la velocidad secundaria en la dirección de la 
superficie del flujo (V,) durante ambos modelos numéricos y pruebas 
experimentales. 
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Nomenclatura 

b es el ancho de la base (m). 

Q es descarga (l/s). 

Lj es longitud de salto hidráulico (m). 

W es componentes de velocidad en dirección horizontal (m/s). 
es componentes de velocidad en dirección vertical (m/s). 

Y) es la viscosidad cinemática (m/s). 

v2 denota la tensión normal de Reynolds (vertical). 

w? denota la tensión normal de Reynolds (horizontal). 

Vr es la velocidad de las corrientes secundarias en dirección 
horizontal (m/s). 

V, es la velocidad de las corrientes secundarias en dirección 
vertical (m/s). 

V es el operador de la nebulosa. 

ij, K son vectores de coordenadas de unidad. 

U es la velocidad en la función de rizo. 

Si, $2 son señales de retorno de dos pulsos adyacentes. 

S1, S2 son parte coherente de la señal. 

Ni, Na son ruido aleatorio. 

m es cotangente de pendiente lateral. 

Fr; es el número de Froude antes del salto hidráulico. 

Y es la profundidad secundaria del salto hidráulico (m). 

SNR es la relación de ruido de la señal. 

Y2/Y, es la relación de profundidad Sequent. 

E,/AE es la pérdida relativa de energía del salto hidráulico. 

RMSE Root-Mean-Square-Error. 

DC es el coeficiente de determinación. 

E es la velocidad del sonido (m/s). 

v es la velocidad del transmisor al receptor (m/s). 

0) es función de flujo. 

Fboppler está cambiando la frecuencia del sonido recibido. 
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Fsourca es la frecuencia de sonido del transmisor. 
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